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ABSTRACT 



We present UVES/VLT observations of the nearby flare star CN Leo covering the Wing-Ford FeH band near 1 jim with high spectral 
resolution. Some of the FeH absorption lines in this band are magnetically sensitive and allow a measurement of the mean magnetic 
flux on CN Leo. Our observations, covering three nights separated by 48 hours each, allow a clear detection of a mean magnetic field 
of Bf x 2.2 kG. The differential flux measurements show a night-to-night variability with extremely high significance. Finally, our 
data strongly suggest magnetic flux variability on time scales as low as 6 hours in line with chromospheric variability. 
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1. Introduction 

Stellar activity can be observed in various tracers from all lay- 
ers of a stellar atmosphere: Doppler imaging techniques reveal 
spots in the photosphere and trace regions of large magnetic flux, 
chromospheric emission shows up in specific lines like e. g. Ha 
or Ca ii H & K originating above photospherically active regions, 
and finally, X-rays trace the hottest material and allow the in- 
vestigation and characterization of coronal properties. In active 
stars, these activity indicators can exceed the levels observed on 
the Sun by several orders of magnitude. 

The strong activity observed in many stars implies strong 
magnetic fields, which are thought to lie at the root of all 
observed activity phenomena. Yet the direct measurement of 
magnetic fields in stellar photospheres is difficult. Most of the 
methods to measure magnetic fields require the detection of 
Zeeman broadening in magnetically sensitive spectral lines. 
Magnetic field strengths up to several kilo-Gauss and large fill- 
ing factors have been measured for a few early M dwarfs and 
T Tauri stars from optical Fei and infra red Tii lines (for an 
overview see iJohns-Rrull & Vaien ti 2000). Spectropolarimetry 
allow s to reconstruct the field geometry for some of these ob- 
jects dDonati et alj |2006|) . However, appropriate atomic lines 
for such measurements either disappear or become more and 
more blended with molecular line s for objects later than spec- 
tral type mid M. iReiner s & Basri (2006) introduced a method 
to measure the mean magnetic flux of these ultracool dwarfs 
from lines of the FeH absorption band around 9950 A. The 
FeH band comprises a large number of absorption lines with 
very different Zeeman sensitivity. Lines with small sensitivity 
can be used to determine the rotation velocity, and magneti- 
cally sensitive lines unambiguously indicate the quantity mag- 
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netic field strength x filling factor (B 0. IReiners & Basril (I2007I) 
successfully measured the mean magnetic field in a sample of 22 
M dwarfs with different activity levels and demonstrated that the 
magnetic flux of their sample objects is strongly correlated with 
Ha emission. 

Stellar activity is inherently variable on time scales from sec- 
onds to decades. Variability on the shortest time scales occurs in 
the form of sudden outbursts or flares, which can be observed 
in several wavelength bands. While this variability has been ex- 
tensively studied through proxy indicators like those mentioned 
above, we tackled the challenge to directly measure magnetic 
field changes and correlate the results with proxies like, e.g., 
Ha emission. For this purpose, we chose to study the flare sta r 
CN Leo (G1406). At a distance of 2.39 pc (iHenrv et aUl2"004T) . 
CN Leo is one of the nearest late-type stars and also among 
the ob jects examined for magnetic fields by IReiners & Basril 
(2007), who estimated an integrated surface magnetic flux of 
Bf « 2.4 kG. The temperature of CN Leo is 2 800 K to 2900 K 
dPavlenko et al . 2006; Fuhrmeister et al. 2005) at a spectral typ e 
of M5.5 or M6.0 dReid et all Il995t iKirkpatrick et alJ fl99lb . 
Further information in particular on the coronal properties of 
CN Leo can be found in Fuhrmeister et al. (2007). CN Leo has 
- fortuitously but somewhat surprisingl y - a rather low apparent 
rotation velocity of v s in i m 3.0km s _1 (Mohantv & Basri 2003; 
IReiners & Basri 2007), thus making it an ideal target for spec- 
troscopic high-resolution studies of stellar magnetic activity. 



2. Observations and data reduction 

CN Leo was observed during three half nights of May 19/20, 
21/22 and May 23/24, 2006 with UVES/VLT. The UVES spec- 
trograph was operated in dichroic mode with spectral coverage 
from 3200 A to 3860 A in the UV and from 6400 A to 10080 A 
with a gap between 8 190 A and 8400 A in the red. We used a 
slit width of 1" resulting in a resolving power of R ~ 40000. 
Here we consider only the far red part of our spectra on the 
upper chip in the red arm. Exposure times varied according to 
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Fig. 1. Variation of the mean magnetic flux Bf over the three 
nights. The relative error on a an individual data point is about 
80 G (see text). 



seeing conditions from 100 s to 300 s, resulting in a total of 
68, 24 and 89 spectra taken in our three half nights. The data 
were processed using the ID L based REDUCE reduction pack- 
age dPiskunov & Valentill2002l) . Since the MIT/LL CCD-20 chip 
shows very strong fringing in the red wavelength region, we per- 
formed a 2D-flatfielding with standard calibration lamp expo- 
sures. The typical signal-to-noise ratio (SNR) of a single spec- 
trum is on the order of 40 around 1 pm. 



3. Method of Analysis 

In order to determine the mean magnetic field strength in 
CN Leo we use the molecular FeH band near 1 pm; this band 
contains a set of narrow, individual absorption lines, some of 
which are magnetically very sensitive, while others are not. 
Unfortunately, so far no Lande g-factors are available for these 
lines. Therefore, magnetic flux measurements are performed by 
comparing FeH lines to spectra of stars with known magnetic 
flux and known activity, i.e., to template stars. The feasibility 
and perfo rmance of this method has been successfully demon- 
strated by IReiners" & Basril (120071) . As a second measurement of 
the magnetic flux, we use the line ratio method suggested by 
iReiners & Basril (|2006) as a consistency check after we calcu- 
lated the magnetic flux with the line-profile fitting method. 

3. 1 . Line profile fitting 

To determine th e mean magneti c field of CN Leo we follow 
the approach by IReiners & Basril d2007l) . who analyzed HIRES 
Keck data of 22 ultracool dwarfs. The basic idea of this method 
consists in modeling the observed spectrum of the target star as 
a linear combination of the spectrum of a magnetically inactive 
star with very low magnetic flux (we chose GJ 1002, M5.5) and 
that of a magnetical ly very active star (we chos e Gl 873, M3.5) 
with Bf * 3.9 kG dJohns-Krull & Valentill2000l) . As actual tem- 
plates we used the Keck HIRES spectra of Gl 873 and GJ 1002. 
A comparison of the UVES and Keck spectra of CN Leo shows 
that the resolving power of the HIRES and UVES spectra are 
very much alike. We therefore employed both data sets without 
any further assumptions on the instrument profiles. 



Specifically, we construct a linear interpolation of the nor- 
malized spectrum ScNLeo °f CN Leo from the normalized tem- 
plate spectra 5 gi 873 and 5gi 1002 through 

ScNLeo - P^G1873 + (1 ~ /?)SgJ 1002- (1) 

The weight p, with < p < 1, yields Bf — p Bf (Gl 873) = 
p 3.9 kG. We fin d the best value of p by x 2 -minimization 
dPress et aT]| 19921) through a comparison of ScNLeo with the ac- 
tually measured (normalized) spectrum. Note that a few steps are 
required before the interpolation of the two template spectra can 
be computed: Both template spectra are adjusted to match the 
target spectrum in line depth and in the overall width of the mag- 
netica lly insensitive lines (for more details, see IReiners & Basril 
2007). The widths of the magnetically insensitive lines are used 
to calibrate the line broadening that affects all spectral lines. This 
global broadening is dominated by stellar rotation and the spec- 
trograph's resolving power. While stellar rotation can obviously 
be assumed constant during all observations, the spectral resolu- 
tion can also be affected by the seeing conditions during observa- 
tion. In a perfect spectrum, for a star with low rotation velocities 
this effect would mimic the effect of extra rotation and should 
show up if we formally allow the fit to also adjust the projected 
rotation velocity v sin / for each spectrum. 

We performed the x 1 -minimization first treating v sin i, Bf 
(i.e., p, see Eq.[T|i, and the line depth a as free parameters. 
However, we did not find any correlation between v sin i and 
the seeing conditions that we measured through the widths of 
telluric lines. We therefore believe that small extra broadening 
due to variable seeing conditions does not significantly affect the 
line widths of the absorption lines, which are dominated by tem- 
perature and stellar rotation. Treating v sin / as a free parameter 
we found v sin / values between 2.0 and 2.5kms _1 , which are 
only marginally different given the comparably low resolution 
(1 resolution element « 7.8kms _1 )Q Since the rotation veloc- 
ity of CN Leo is not variable on the timescales of our measure- 
ments, we took the mean value of all values of v sin / from the 
fit, v sin / = 2.14km s _1 (incidentally, the value of the mean co- 
incides with the value of the median, so it is not necessary to 
decide which measure is best here). We then repeated the fit, 
now with a fixed value of v sin / = 2.14km s . This procedure 
did not significantly change the derived results when compared 
to the results achieved with v sin / as a free parameter. In par- 
ticular, the derived Bf values do not vary significantly, and we 
proceed using the results achieved with the fixed value of v sin i, 

A total of 1 8 1 spectra of CN Leo were recorded during our 
campaign and were used to derive mean magnetic fields as de- 
scribed above. In Fig.Q] we plot all our Bf measurements versus 
time. The typical uncertainty of our Bf measurements is approx- 
imately ABf = 80 G; this value represents the 1-cr uncertainty 
from the y 2 -fit, i.e., the values of Bf for which y 2 < v 2 . +1 

J ^1 — l A min 

dPress et al.ll 19921) . We emphasize that this error is of purely sta- 
tistical nature and cannot be identified as the true uncertainty of 
the absolute measurement of the magnetic flux Bf. For the rela- 
tive measurement, however, this is the uncertainty that we adopt. 
We note that our fitting procedure started with an initial guess 
of Bf = 2500 G for every spectrum, and that the algorithm did 
significantly modify this initial guess. We also tried the effect 
of a different but fixed value of v sin / and found that a different 
v sin i slightly modifies the absolute values of all measurements, 
but does not change our results on the differential measurements. 

1 Note that v sin i is the broadening required to match the template 
spectra and the spectra of CN Leo, thus it is only a lower limit for the 
true projected rotation velocity of CNLeo. 
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3.2. Analysis of spectral differences 



Table 1. Mean magnetic flux, standard deviation, error of the 
nightly mean, number of observations per night for all three 
nights; the other columns contain results of non-parametric sta- 
tistical tests (see text). 



Night 


Bf± 


c [G] 


N 


KS1 


KS2 


RSI 


RS2 


1 


2148 ± 


45 ± 5.5 


68 










2 


2248 ± 


47 ± 9.6 


24 


7 kt 11 









3 


2275 ± 


48 ±5.1 


89 


3 10- 28 


8.4 10" 3 





7.7 10-' 



Is there any true variation of Bf visible in our data? The dif- 
ference in the magnetic flux between the individual spectra with 
the maximum and the minimum flux is ss 350 G, i.e., less than 
10 % of the difference in flux between the two template stars of 
« 3900 G. At the spectral type of CN Leo, the line depth of the 
magnetically most sensitive lines varies by » 15 % of the con- 
tinuum value (e.g. t he line at 9948 A varie s between ~ 0.58 and 
~ 0.74, cp. Fig. 6 in lReiners & Basril 120071) . Thus, in order to de- 
tect a magnetic flux variation of 350 G, a variation on the order 
of 1.5 % has to be detected in this line, i.e., a SNR on the order 
of 70 is required. Because of the short individual exposures such 
a SNR is not available in our spectra, and therefore a direct com- 
parison cannot be expected to yield a significant difference in 
the Zeeman signatures directly visible in the spectra. However, 
our x 1 -At uses simultaneously the information contained in many 
lines and hence is superior to a direct comparison of line depths. 

Nonetheless, our^ 2 -fits could be fooled by a variety of ef- 
fects that are difficult to control if the results cannot be inde- 
pendently verified. In order to achieve such a consistency check, 
we constructed two mean spectra by averaging all spectra from 
the first night (68) and all spectra from the third night (89). The 
mean values of Bf for the three different nights are given in 
TableQ] Between the first and the third night ~Bf differs by 125 G. 
The summation of all spectra in the first and third night with a 
SNR of typically ~ 40 each yields average spectra with a SNR on 
the order of ~ 300. Since a SNR on the order of 200 is required 
to detect a difference of 125 G in the 9948 A lines between the 
two mean spectra (see above), the difference in the magnetically 
most sensitive lines should become visible in such a comparison. 

In Fig.|2l we show the mean spectra for the first and third 
nights. The magnetically very sensitive line at 9948 A is marked 
with a "Y", two relatively insensitive lines are marked with an 
"N". As expected from our estimate, a small difference becomes 
visible in the magnetically sensitive line, while insensitive lines 
do not show systematic variation. 



3.3. Magnetically sensitive line ratios 

As a second consistency check we use the two mean spectra for 
the first and thi rd nights and calculate the two line depth ratios 
as suggested bv lReiners & Basril ([2006). Each ratio is calculated 
from measuring the depth of two neighbored lines, one of which 
is Zeeman sensitive and one is insensitive. We us e the two ratios 
9949. 1 1/9948.13 and 9956.77/9958.25 as in iReiners & Basril 
(2006), both ratios should be larger at high magnetic flux. The 
re sults are given in Tabl e[2] and should be compared to Fig. 7 
of I einers & Basril ((2006). That figure shows that indeed the 
9949.11/9948.13 ratio is more sensitive to magnetic flux vari- 




9946 9948 
Wovelength[A] 



Fig. 2. Co-added spectra from the 68 spectra of the first night 
(red) and from the 89 spectra of the third night (black). At the 
achieved SNR « 300 Zeeman broadening becomes visible in the 
most sensitive lines. Y: magnetically sensitive line, N: magneti- 
cally insensitive lines. 



ation than the 9956.77/9958.25 ratio is. Both ratios and their 
variation are in good agreement with our results from ^-fitting. 



Table 2. Line depth ratios as introduced in IReiners & Fjasri 
(120061) . Both ratios are consistent with our results (see text). 





9949.11/9948.13 


9956.77/9958.25 


S/Nightl 


1.48 ± 0.01 


1.44 ±0.01 


S/Night3 


1.59 ±0.01 


1.45 ±0.01 



4. Results, Discussion and Conclusions 

As a first result of our analysis of 1 8 1 spectra of CN Leo we note 
that magnetic flux is consistently detected in each UVES spec- 
trum. At first sight the data are consistent with the assumption of 
a constant magnetic flux during each night, while the magnetic 
flux does appear to vary from night to night. In TableQ] we give 
the derived mean values of B f and their standard deviation for 
all three nights. The standard deviation is consistent with the es- 
timated uncertainty of an individual observation; in fact, it is a 
even little smaller, which could mean that the differential mea- 
surements are somewhat more precise than estimated. 

The uncertainty in the absolute value of the field is domi- 
nated by our ignorance of the splitting patterns of the magneti- 
cally sensitive FeH lines. Our method only employs a linear in- 
terpolation between the two spectra of the template stars, i.e., it 
assumes that the splitting patterns are the same (B is constant), 
and that only the filling factor / is variable. This is probably not 
the case, although the change in filling factor is usually assumed 
to be the dominant factor. At any rate, we can not yet calculate 
differe nt splitting patterns and thus, following IReiners & Basril 
(120071) . we estimate the uncertainty on the absolute value to be 
on the order of a few hundred Gauss. However, here we are in- 
terested in the accuracy of the absolute measurements only to 
a lesser extent, rather are we interested in the differential mea- 
surements comparing the values of Bf during several nights and 
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Fig. 3. Magnetic flux (left panels) and Ha emission (right panels) 
for first night (lower panels) and third night (upper panels); in 
each panel, data for the respective (half-)nights are shown for 
the first half (solid histograms) and second half (dashed). Note 
the correlated change in magnetic flux and Ha emission during 
the first and third half night. Eight data points in a large flare 
during the first half of the first night are not shown because of 
saturation of the Ha line core. 



within a given night. Since the magnetically insensitive lines stay 
constant, they provide a very accurate calibration, particularly 
when intercomparing the exposures. Thus, we estimate that our 
uncertainties in the individual differential measurements are in- 
deed on the order of the uncertainty of the j^-fit, i.e., 80 G. 

Inspecting the nightly mean magnetic field values and their 
dispersion (TableQ}, the difference in these means appears to be 
significant only on the 1-2 cr level when compared to the vari- 
ance of the distribution. We checked the nature of the distribu- 
tion of our nightly measurements around the respective nightly 
means and found very good agreement with a Gaussian distri- 
bution and therefore the accuracy of the mean ought to be much 
smaller than the width of the observed distribution. Computing 
the uncertainty of the nightly means, we find values reduced by 
y/N (cf., TabQ]l and the differences between the individual be- 
come very significant indeed. We verified the existence of night- 
to-night magnetic field variations also by analyzing the actu- 
ally measured distributions of Bf using two non-parametric test 
methods. We specifically applied the Wilcoxon rank sum test 
(denoted by RS in Table [TJ and the two sample Kolmogorov- 
Smirnov test (KS) to investigate whether the measured Bf dur- 
ing our three observing nights are drawn from the same parent 
population (KS 1 and RS 1 indicate comparison to the first night, 
KS2 and RS2 to the second). The derived probabilities from 
the KS- and RS-tests that the null hypothesis "no night-to-night 
magnetic flux change" is correct are also listed in Table Q] As is 
clear from TableQ] the field measurements for the nights 1 and 2 
and 1 and 3 differ in an extremely significant fashion, and even 
the measurements between nights 2 and 3 appear to differ at a 
confidence level in excess of 0.99 independent of the test statis- 
tics used. In summary, the above statistical analysis strongly sug- 
gests magnetic flux changes on time scales of 48 hours. 

How about the presence of flux changes on shorter time 
scales? In order to address this issue we subdivided the avail- 
able magnetic flux measurements for each of our three nights 
into two halves and applied the same tests performed for night- 



to-night variability to assess any intra-night variability. While we 
cannot claim any intra-night variability during night 2 (with only 
24 Bf measurements available), we find significant intra-night 
variability for night 1 and night 3. The variability is such that the 
magnetic flux during night 1 appears to be increasing, while it is 
decreasing during night 3. Interestingly, the intra-night variabil- 
ity in chromospheric emission lines follows the same pattern as 
the magnetic flux. We measured the c hromospheric flux fro m the 
Ha lines contained in our spectra (see lReiners & Basril2007f) and 
calculated L Ha /L\, ] with the same value of Lboi f° r a ll spectra. 
The change in magnetic flux and Ha emission is directly visi- 
ble in the distribution of magnetic flux and Ha emission during 
the first and second quarters of nights 1 and 3 (first and second 
half of our observation in each night). The distribution for nights 
1 and 3 are shown in the lower and upper panels of Fig. [3] re- 
spectively. We find the difference between the two distributions 
shown in each panel of Fig.[3]is statistically highly significant. 

These trends are verified by formally describing the nightly 
magnetic flux measurements as a linear function with time. 
Linear regression analysis shows that including a linear trend 
significantly better describes the data as opposed to the assump- 
tion of a constant magnetic flux; the values of A% 2 are 6.0 and 
7.3 for nights 1 and 3 respectively. We find a positive flux change 
during night 1 and a negative change during night 3; the ob- 
served flux changes to an amount of « 50 G over a time scale of 
6 hours. While we cannot entirely exclude artifacts of the fit pro- 
cedure as the cause for the intra-night variations, two facts make 
us confident that the observed changes are of a physical origin: 
First, the observed changes are of opposite sign during night 1 
and 3. Second, the magnetic flux changes during those nights 
correlate with corresponding activity changes during the same 
nights, which we have verified in the Ha lines contained in our 
spectra (see Figf3]). We note that our spectra contain a multitude 
of chromospheric emission lines, which are correlated with Ha 
em ission. These will be an alyzed in a subsequent publication. 

Rein ers & Basril (|2007) measured Bf = 2.4 kG in a spec- 
trum of CN Leo taken in 2005. This result is on the upper end 
of the values we found during our observations. However, their 
spectra were not taken at the same instrument, which introduces 
uncertainties due to the different calibrations. Hence their result 
is consistent with our measurements. 

In summary, we have successfully measured a mean mag- 
netic field of Bf x 2.2 kG on the flare stars CN Leo using the 
FeH bands near 1 /mi. Our monitoring campaign shows definite 
magnetic flux changes of ~ 100 G on a time scale of 48 hours. 
Further, the data of the individual (half) nights suggest flux 
changes on even shorter time scales, which are correlated with 
the general activity level of CN Leo as measured in Ha emis- 
sion. Further monitoring of magnetic fields and their variations 
in CN Leo and similar stars appears very promising. 
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